Volume fluctuations and higher order cumulants of the net baryon number 
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We consider the effect of volume fluctuations on cumulants of the net baryon number. Based 
on a general formalism, we derive universal expressions for the net baryon number cumulants in 
the presence of volume fluctuations with an arbitrary probability distribution. The relevance of 
these fluctuations for the baryon-number cumulants and in particular for the ratios of cumulants is 
assessed in the Polyakov loop extended quark-meson model within the functional renormalization 
group. We show that the baryon number cumulants are generally enhanced by volume fluctuations 
and that the critical behavior of higher order cumulants may be modified significantly. 



I. INTRODUCTION 

One of the goals of the experiments with ultrarelativis- 
tic heavy ion collisions at SPS, RHIC and LHC energies is 
to probe the phase structure of strongly interacting mat- 
ter and, in particular, to identify the deconfinement and 
chiral restoration transitions. In this context, the fluc- 
tuations of conserved charges may serve as a pertinent 
probe. 

Fluctuations of the net baryon number and electric 
charge may provide an experimental signature for the 
hypothetical chiral critical endpoint [TJ [5] . Moreover, as 
recently noted j3HH] , such fluctuations are also of inter- 
est at small baryon densities, since they reflect the critical 
dynamics of the underlying 0(4) transition, expected in 
QCD in the limit of massless light quarks [T21 [T3]. In- 
deed, it was demonstrated that higher order cumulants 
change sign in the crossover region of the QCD phase dia- 
gram jHIHllS]- Thus, the observation of a strong suppres- 
sion the higher order cumulants may be used to identify 
the chiral crossover transition in experiment. 

The first measurements of fluctuations of the net 
baryon number, more precisely of the net proton num- 
ber 1 , in heavy ion collisions at RHIC were obtained by 
the STAR Collaboration [T5]. The analysis of cumulants 
of the fluctuations and of the probability distributions 
confirmed, that the hadron resonance gas (HRG) model, 
which yields a quantitative description of particle yields 
in heavy ion collisions |16) , provides a useful reference for 
the non-critical background contribution to the charge 
fluctuations j4,. Thus, critical fluctuations related to the 
dynamics of the chiral transition should be reflected in 
deviations of the measured net charge fluctuations from 
the HRG baseline. In this context, higher order cumu- 
lants are of particular interest [H [10] . 

A detailed analysis of experimental data on moments 
of net proton number fluctuations and their probability 
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1 As shown recently in Kef. 1141 , the net nucleon number cumu- 
lants can, to a good approximation, be deduced from the mea- 
sured net proton cumulants. 



distributions indeed exhibit deviations from the HRG. 
To verify the origin of these deviations, one must iden- 
tify and assess effects, unrelated to the critical dynamics, 
which can influence the charge fluctuations. For instance, 
it was recently argued that constraints, owing to the con- 
servation of the total baryon number in nucleus-nucleus 
collisions, modify the non-critical background contribu- 
tions to higher order cumulants of the net proton number 
fluctuations [17] , 

In this paper we study volume fluctuations as a fur- 
ther possible source of non-critical fluctuations, not ac- 
counted for in the HRG model results. We first present 
a transparent derivation of the cumulants of net baryon 
number, including the effect of volume fluctuations. The 
resulting cumulants are expressed in terms of cumulants 
of the net baryon number distribution at fixed volume 
and cumulants of the probability distribution for volume 
fluctuations. We also provide a more formal derivation, 
making use of the cumulant generating functions. We 
stress that the final expressions are general, independent 
of the probability distributions for net baryon number 
and volume. The only assumption made is that the two 
sources of fluctuations are independent. 

Furthermore, we assess the effect of volume fluctua- 
tions on the kurtosis and on ratios of higher order cumu- 
lants of the net baryon number in the vicinity of the chi- 
ral transition within the Polyakov loop extended quark- 
meson model. Here we employ the functional renormal- 
ization group to properly account for critical fluctuations. 

The paper is organized as follows: In the next section 
we obtain the corrections due to volume fluctuations to 
the first four moments of the net baryon number fluctu- 
ations. In Section III we derive a general expression for 
the corrected cumulants, valid to any order, obtained us- 
ing the cumulant generating functions. In Section IV we 
illustrate the role of volume fluctuations with a numerical 
study and finally in Section V we state our conclusions. 



II. HEURISTIC APPROACH 

Consider a fixed volume V, where the net baryon 
number B fluctuates with the probability distribution 
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P(B,V). The n-th order moments of the net baryon 
number are then defined by 



(B n )v = 



OO 

£ 



B n P{B,V). 



(1) 



It is convenient to introduce reduced cumulants, corre- 
sponding to the net baryon number fluctuations per unit 
volume. The first four reduced cumulants are 
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K4(T,fl) = 
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where 5B = B — B. The cumulants to leading 

order, independent of the volume V. In the following we 
neglect subleading surface effects, which could lead to a 
residual volume dependence of the cumulants. 

The volume dependence of the moments follows from 
([2| and reads 

(B) v = Kl V, 
(B 2 ) v = k 2 V + k\ V 2 , 

(B 3 ) v = k 3 V + 3k 2 k 1 V 2 + k 3 V 3 , (3) 
(B A ) V = n 4 V + (4k 3 ki + 3k 2 .) V 2 + Qk 2 k\ V 3 + k\ V 4 . 

The coefficients in Eq. ([3| are those of the Bell polyno- 
mials. 

As an illustrative example, we consider the hadron res- 
onance gas. In this model, the net baryon number fluc- 
tuations are given by the Skellam distribution E] and 
the corresponding cumulants are particularly simple 2 

1 ,~ - , (HRG) 1 



(HRG) 



2n+i -y{B x -B^), *£™> = V (B 1+ B^), (4) 



where B\ = [B\) is the mean number of baryons and 
B-x = (B-i) that of anti-baryons in V. The correspond- 
ing moments are obtained by inserting the cumulants Q 
in ((§. 

We now allow for fluctuations of the volume. To this 
end, we introduce the volume probability distribution 
■p(V), the corresponding moments 



(V n ) = / V n P(V)dV, 



(5) 



and the reduced cumulants of the volume fluctuations, 
v n . The latter arc defined as in Eq. ^ with the replace- 
ments V — > (V) and B — > V. Thus, e.g. v% = 1 and 

v 2 = ((V 2 )-(V) 2 )/(V). 



In the presence of volume fluctuations the moments of 
the net baryon number are given by 

r. OO 

(B n )= dVV(V) B n P(B,V) 

= J dVV(V)(B n ) v . (6) 



It is now straightforward to compute the reduced cumu- 
lants, including the effect of volume fluctuations. Using 
Eqs. {2), (H and ©, we find the general relations 
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K 2 + K x V2i 

n 3 + 3k 2 kxv 2 + k\v 3 , (7) 
k 4 + (4k 3 ki + in\)v 2 + 6k 2 k 2 v 3 + nfvi, 



which are valid for arbitrary probability distributions, 
provided the fluctuations in baryon number and volume 
are independent. We note that the form of Q is deter- 
mined by the volume dependence of the moments d3]). 
Hence, the coefficients in (JtJ) are also given by the Bell 
polynomials. 



III. GENERAL DERIVATION 

In the previous section, we explored the effect of vol- 
ume fluctuations on the fluctuations of the net baryon 
number for the first few cumulants, where explicit calcu- 
lations are tractable. In the following we derive a general 
expression for the cumulants, under the assumption that 
the fluctuations of baryon number and volume are inde- 
pendent. 



A. Formalism 

In general, the probability distributions introduced in 
section II are characterized by the corresponding cumu- 
lant generating functions 3 , 

OO 

X B (t) = hi P(B)exp(Bt), (8) 

B=-oo 

OO 

X V (s) = In J dVV(V)exp(Vs). (9) 



The cumulants are obtained by expanding \ B and \ v in 
series about the origin. The additivity of cumulants and 



2 The normalization of the generalized susceptibilities given in [6] 
differs from the cumulants used here by a factor T 3 . 



3 We assume that the integrals in and |9| converge for t and s 
in an interval around the origin, so that the cumulant generat- 
ing functions exist |18) . Moreover, we assume that the relevant 
cumulants exist. 



3 



thermodynamic principles imply that 

x B (t) = v-( B (t), (io) 

where ( B is a volume-independent function. In fact, ( B 
is the generating function for the reduced cumulants, de- 
fined in Eqs. 

d n 



(11) 



Similarly, we find for the reduced cumulants of volume 
fluctuations 
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(V) dtf 



;X V (t) 



(12) 
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Our aim is to compute cumulants of the net baryon 
number including the effects of volume fluctuations. 
These cumulants are obtained from the cumulant gen- 
crating function 

B (t) = In J dVV(V) P(B, V)e Bt . (13) 



Using Eq. ([10]) wc find 

J2P(B,V)e Bt =e v < B ®, (14) 

B 

and consequently 

cj> B {t) = ]n f dVV{V)e v ^ B{t) . (15) 



A comparison with the definition of the cumulant gen- 
erating function (|9j), yields 

^(*)=X V [C B (*)]. (16) 

This is the general form of the cumulant generating func- 
tion for fluctuations of the net baryon number, including 
the effect of volume fluctuations. The corresponding re- 
duced cumulants are given by a Taylor expansion of (f> B (t) 
about t = 0, 



'(*) 



(17) 



1 d n 

Cn ~Jv)d¥' 

We note that since ( B (t = 0) = 0, no further normaliza- 
tion is needed in the calculation of the cumulants. 

Using Faa di Bruno's formula [19], we obtain a closed 
form expression for the cumulants, 



v n— z-t-1 



), 



(18) 



where B n ^ are Bell polynomials. This equation confirms 
and extends our previous results for the first four cumu- 
lants, given in Eq. Q. Thus, for an arbitrary probabil- 
ity distribution for the fluctuations of net baryon num- 
ber as well as for the fluctuations of the volume, Eq. ( 18 ) 



yields cumulants that can be confronted with experiment. 
Conversely, given a model for the volume fluctuations, 
Eq. ( 18 1 can be used to extract cumulants of the net 



B. Gaussian volume fluctuations 

If the probability distribution for volume fluctuations 
is approximately Gaussian, i.e. v n -C 1 for n > 2, the 
higher order cumulants are simplified considerably 



= K 3 + 3kiK 2 U2, 



c 4 



(3k 2 + 4ki« 3 ) V2, 

(IOK3 + 15K2K4 + 6K1K5) 1>2, 



«4 H 

(35k? + 56K3K5 + 28^2^6 + 8K1K7) V2- 



(19) 



In the section IV we use this form to explore the effect 
of volume fluctuations on the cumulants of net baryon 
number. 



IV. NUMERICAL RESULTS IN THE PQM 
MODEL 



We illustrate the influence of volume fluctuations on 
net baryon number fluctuations, within a model calcula- 
tion. Of particular interest is the modification of higher 
order moments near the chiral crossover transition. We 
adopt the Polyakov loop-extended Quark Meson Model 
(PQM) and compute the cumulants in a non-perturbative 
scheme, the functional renormalisation group. Details on 
the calculations and on the derivation of the net baryon 
number fluctuations can be found in Ref. [Uj. In this 
exploratory calculation, we consider only Gaussian fluc- 
tuations of the volume for a system at zero baryon chem- 
ical potential. In this case, all odd cumulants of the net 
baryon number vanish, i.e. c 2 „ + i = K2n+i = 0. 

In Ref. [8] it was shown that near the chiral crossover 
transition, higher cumulants of the net baryon number 
(n > 4) differ considerably from the predictions of the 
HRG model. In particular, it was suggested that nega- 
tive values of kq and ks could be used to map out the 
chiral phase boundary. This potential signal for the QCD 
phase transition may be affected by volume fluctuations. 



Indeed, the second term in Eq. ( 19 1 yields a positive con- 
tribution to cq. The strength of this contribution is di- 
rectly proportional to the second cumulant i; 2 . 

In general, volume fluctuations are difficult to assess. 
In heavy ion collisions they depend on the centrality of 
the collision, on the definition used to fix the number of 
participants and on the kinematical window, where the 
fluctuations are measured. Thus, u 2 is specific to a given 
experimental setup. We choose t>2 = 0.6, employed in 
Ref. [2D]. 

It is useful to consider ratios of cumulants, 



Rr, 



Cn 



(20) 



baryon number in a fixed volume. 



where the dependence on the (mean) volume of the 
subsystem cancels. In Figs. [T] and [2] we show the 
temperature-dependence of the i?4^ and i?6,2 obtained 
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FIG. 1: The ratio 7?4,2 as a function of temperature com- 
puted in the PQM model at vanishing chemical potential with 
and without volume fluctuations. The probability distribu- 
tion for volume fluctuations is assumed to be Gaussian with 
variance V2 = 0.6, see the text for details. 



in the PQM model for finite and vanishing vi- Moreover, 
in Fig. [3j the dependence of the ratios R n ^ m on v 2 at the 
crossover transition temperature is illustrated. 

These results show that volume fluctuations tend to 
suppress the signature of the chiral transition in the cu- 
mulants of net baryon number. Here, the ratio i?6,2 seems 
to be particularly sensitive. Consequently, the usefulness 
of fluctuations of conserved charges as a probe of crit- 
icality in heavy ion collisions, depends crucially on the 
possibility to control volume fluctuations. 

V. CONCLUSIONS 

We have studied the influence of volume fluctuations 
on the properties of cumulants of net charge distributions 
in heavy ion collisions. In particular, we have computed 
the contribution of volume fluctuations to cumulants of 
the net baryon number. 

In a heuristic approach we showed explicitly how the 
corrections due to volume fluctuations arise. The ex- 
pressions, which hold for arbitrary probability distribu- 
tions, were confirmed and extended in a general formal- 
ism, where we employed the cumulant generating func- 
tions to obtain a closed form for the cumulants, including 
volume fluctuations. 

We assessed the contribution of volume fluctuations for 
the kurtosis as well as for ratios involving higher or- 
der cumulants, viz. i?6,2 and i?8.2, in the Polyakov loop 
extended quark-meson model. A non-perturbative treat- 
ment of fluctuations was obtained by employing the func- 
tional renormalisation group. We focused on the struc- 
ture of the ratios of cumulants in the vicinity to the chiral 
crossover transition, assuming that the probability distri- 
bution of volume fluctuations is approximately Gaussian. 




0.7 0.75 0.8 0.85 0.9 0.95 1 1.05 
T/T pc 

FIG. 2: Same as Fig. [lj but for i? 6 , 2 . 

Finally, we showed that phenomenologically relevant 
ratios of cumulants of the net baryon number are en- 
hanced by volume fluctuations. Consequently, for a suffi- 
ciently large V2, the negative structure of these ratios, ex- 
pected in heavy ion collisions owing to the critical chiral 
dynamics, may be modified significantly by volume fluc- 
tuations. It follows that fluctuations of conserved charges 
in heavy ion collisions can provide a robust probe of the 
chiral phase boundary, if good control of volume fluctu- 
ations can be achieved. 
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FIG. 3: The ratios R„ t m at the chiral crossover temperature 
T pc , obtained in the PQM model, as functions of the second 
cumulant of the volume fluctuations, 1)%. 
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